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ABSTRACT 



In general there are a large number of light scalar fields in the theories going 
beyond standard model, such as string theory, and some of them can be taken 
q_) ■ as the candidates of curvatons. For simplicity, we assume all of curvatons 

have the same decay rate and suddenly decay into radiation at the same time. 
In order to distinguish this scenario from the more general case, we call it 
"N-vaton" . We use 5J\f formalism to calculate the primordial power spectrum 

t*-. \ and bispectrum in N-vaton model and investigate various bounds on the non- 

Gaussianity parameter f^L- A red tilted primordial power spectrum and a 
large value of /jvl can be naturally obtained if the curvature perturbation 

(-~. generated by inflaton also makes a significant contribution to the primordial 

power spectrum. As a realistic N-vaton model, we suppose that the axions in 
the KKLT compactifications of Type IIB string theory are taken as curvatons 
and a rich phenomenology is obtained. 
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1 Introduction 

Most inflation models predict a nearly Gaussian distribution of the primordial curvature 
perturbation. Deviations from an exactly Gaussian distribution are characterized by a 
dimensionless parameter f NL [1]. In the case of single field inflation model j'nl ~ 0{n s — 1) 
[2], which is constrained by WMAP (n s = 0.960^013) [3] to be much less than unity. A 
Gaussian distribution of the primordial curvature perturbation is still consistent with 
WMAP five-year data [3]: 

- 9 < f^t < 111 and - 151 < f e N q f < 253 (95%CL), (1.1) 

where "local" and "equil" denote the shapes of the non-Gaussianity. In [4] the authors 
reported that a positive large non-Gaussianity 

27 < f^f < 147 (1.2) 

is detected at 95% C.L.. Planck is expected to bring the uncertainty of /^£ a/ to be less than 
5 [5]. If a large value of Jnl is confirmed by the forthcoming cosmological observations, 
the simplest model of inflation is ruled out and some very important new physics of the 
early Universe will be showed up. 

In general a large number of light scalar fields are expected in the theories beyond the 
standard model, such as string theory. The consistent perturbative superstring theory 
can only live in ten-dimensional spacetime. To connect string theory with experiments, 
string theory must be compactified on some six-dimensional manifold and many dynam- 
ical moduli fields emerges in four dimensions. The typical number of moduli fields is 
N ~ (9(10 2 — 10 3 ). One can expect that the expectation values of some of these scalar 
fields are displaced from the minimum of their potential due to the quantum fluctuations 
during inflation. Usually these scalar fields are subdominant during inflation and their 
fluctuations are initially of isocurvature type. After the end of inflation they are supposed 
to completely decay into thermalized radiation before primordial nucleosynthesis and thus 
the isocurvature perturbations generated by them are converted to be a final adiabatic 
perturbations. These scalar fields are called curvatons. 

The curvaton mechanism to generate an initially adiabatic perturbation deep in the 
radiation era is proposed in [6-9]. The primordial curvature perturbation in curvaton 
model with single curvaton has been discussed in [6-13]. If many curvatons make contribu- 
tions to the primordial density perturbation, the calculation becomes much more compli- 
cated [14-16]. Since curvaton model can give a large positive local- type non-Gaussianity, 
recently some topics related to curvaton model are discussed in [17-24]. 
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In single-curvaton model Jnl is inverse proportional to the fraction of curvaton energy 
density in the energy budget at the epoch of curvaton decay. Smaller the energy density of 
curvaton, larger the non-Gaussianity. Since the curvaton mass is smaller than the Hubble 
parameter H* during inflation, or equivalently its Compton wavelength is large compared 
to the curvature radius of the de Sitter space H~ , the gravitational effects play a crucial 
role on the behavior of curvaton field in such a scenario. The typical energy density of 
curvaton field in such a background is roughly H%, which leads to an upper bound on 
/nl [17]: /ml < 522 -r* (up to an order one coefficient), where r is the tensor-scalar ratio. 

Since multiplicity of scalar fields is generally expected, we focus on the multi-curvaton 
scenario in this paper. Here we consider a special case in which all of curvatons have the 
same decay rate and their masses are larger than the decay rate. In order to simplify the 
calculation of the primordial curvature perturbation, we assume all of curvatons suddenly 
decay into radiation at the same time. We give a name, "N-vaton", to this scenario. As 
a realistic N-vaton model, the axions in the KKLT compactification of Type IIB string 
theory are suggested to be curvatons. Based on the random matrix theory, the mass 
spectrum of axion obeys the Marcenko-Pastur law and a rich phenomenology of this 
model is shown up. 

Our paper is organized as follows. In Sec. 2, we use the SJ\f formalism [25-27] to 
calculate the primordial curvature perturbation on large scales in N-vaton model. The 
various bounds on f^L are discussed in Sec. 3. In Sec. 4 we consider a more general case 
where the curvature perturbation generated by inflaton cannot be ignored and we find 
that the spectral index of primordial power spectrum can be red-tilted naturally. In Sec. 
5, we propose a realistic N-vaton model in which the axions in the KKLT compactification 
of Type IIB string theory are taken as curvatons. At the end we give some discussions on 
N-vaton model in Sec. 6. 

2 Primordial curvature perturbation 

In this paper we consider that inflaton <p and curvatons cr, are decoupled to each other. 
The action takes the form 
Ml 



S = — P - J d^Xy^JR + f d 4 Xyf^) 



z j=l z 



(2.1) 



where M p = 2.438 x 10 18 GeV is the reduced Planck scale and the potential V(4>, <7j) is 

given by 

1 N 
^(0,a,) = V(0) + -]Tm^ 2 . (2.2) 

i=l 



During inflation the total energy density is dominated by inflaton potential V(<f>) and 
the dynamics of the system is described by the equation of motion of inflaton and the 
Friedmann equation: 

+ 3^0+^ = 0, (2.3) 

We also define some slow-roll parameters, such as 

2 V V(<j>) J ' * V(0) ( ' 

If e <C 1 and \rj\ <C 3, inflaton slowly rolls down its potential. 

In this paper we expand any field or perturbation at each order (n) as follows 

OO 1 

C(t,x) = C (1) (*,x) + £-rC (n) (*,x). (2.6) 

n=2 n - 

We assume that the first-order term (^ is Gaussian and higher-order terms describe the 
non-Gaussianity of the full nonlinear (. Working in the framework of Fourier transforma- 
tion of C, the primordial power spectrum Vq is defined by 

(C(k 1 )C(k 2 )) = (27r) 3 V c (h)6 3 (k 1 + k 2 ), (2.7) 

and the primordial bispectrum takes the form 

(C(k 1 )C(k 2 )C(k 3 )) = (27r) 3 J B c (k 1 , k 2 )5 3 (k! + k 2 + k 3 ). (2.8) 

The amplitude of the bispectrum relative to the power spectrum is parameterized by the 
non-Gaussianity parameter f NL , i.e. 

B c {kt, k 2 ) = ^MVcik^Vcih) + 2 perms]. (2.9) 

The primordial density perturbation can be described in terms of the nonlinear cur- 
vature perturbation on uniform density hypersurfaces [28] 

C(t,x) = M/-(;,x) + -/ -^, (2.10) 

3 Jp(t) p + P 

where Af — J Hdt is the integrated local expansion, p is the homogeneous density in the 
background model, p is the local density and p is the local pressure. 



After inflation inflaton decays into radiations which dominate the total energy density 
of our universe. In the radiation dominated era the Hubble parameter H goes like ~ 
a~ 2 . Once the Hubble parameter drops below the mass of curvaton field the field starts 
to oscillate. Nonlinear evolution of the values of curvatons on large scale is possible if 
the potential of curvatons deviates from a purely quadratic potential away from their 
minimums [29,30]. Thus, in general, the initial amplitude of curvaton oscillations <Tj. is 
some function of the field value o^* at the Hubble exit lj: 

0i,o = 9i{vi,*)- (2.11) 

When the curvaton starts to oscillate about the minimum of its potential, but before 
it decays, it behaves like pressureless dust (p CTi o ~ a -3 ) and the nonlinear curvature 
perturbation on uniform-curvaton density hypersurfaces is given by 

C^(t,x) = ^(t,x) + p° (t ' X) ^. (2.12) 

The curvaton density on spatially flat hypersurfaces is 

p.JsM=o = e K ^p^ . (2.13) 

The quantum fluctuations in a weakly coupled field, such as curvaton, at Hubble exit 
during inflation are expected to be well described by a Gaussian random field [31]. So we 
have 

0i,* = o\* + <foi,*, (2-14) 

without higher-order nonlinear terms. During the curvaton oscillation we expand the 
energy density p ffio = \mfal and ( aio to second order: 



P<?i, o P&i,o 



l + 2X i +(l + h i )Xi 



(2.15) 



2 
where p aio = |m?ff? , a i;0 = ^(ct^), and 



c^ = c£ + k£, (2.i6) 



6a W 
X % = -^, (2.17) 

0"t,o 

hi = ^#. (2.18) 



x In this paper the subscript * denotes the quantity evaluated at the Hubble exit. 



Here prime denotes the derivative with respect to 01*. Order by order, from Eg. (12. 131) we 
have 

&l = \*i> (2-19) 

C& = -f(l-fc)($„) a . (2-20) 

In N-vaton model, we assume that the curvatons have the same decay rate IV . When 
the Hubble parameter drops below IV, all of curvatons decay into radiations. In order to 
get analytic expressions, we work in the sudden-decay approximation which means that all 
of curvatons suddenly decay into radiations at the time tjj when H — Y ff . For simplicity, 
we assume m s > IV for i — 1,2, ...,N and then all of curvatons begin oscillating before 
they decay. 

The curvatons-decay hypersurface is a uniform-density hypersurface and thus from 
Eq. (12.10p the perturbed expansion on this hypersurface is 5Af = (, where ( is the total 
curvature perturbation at curvatons-decay hypersurface. Before the curvatons decay, 
there have been radiations produced by decay of inflaton. Since the equation of state of 
radiation is p r = p r /3, the curvature perturbation related to radiations is 

( r = ( + ]\n^. (2.21) 

4 p r 

The curvatons behave like pressureless dust (p ai = 0) and thus 

C„ = C + ^ln^. (2.22) 

In the absence of interations, the curvature perturbations ( r and ( (7i o are conserved re- 
spectively and the above two equations can be written as 

p r = p r e 4 ^\ (2.23) 

An, D = P^ o e 3iC °"°-°. (2.24) 

At the time of curvatons decay, the total energy density ptot is conserved, i.e. 

N 

p r (t D ,x) +^p CTio (t D ,x) = ptot(t D ). (2.25) 

Requiring that the total energy density is uniform on the decay surface, we have 

N 

(1 - n fftD ) e 4 ^-« + J2 ^De 3{( ^-° = 1, (2.26) 



where fi ai ,D = Pa iD l ' ptot is the fraction of curvaton energy density in the energy budget 
at the time of curvaton decay, and 



N 



Qff,D = 7 , Qcr,.,D- 



(2.27) 



i=l 



Actually ( r is generated by the fluctuation of inflaton during inflation, namely ( r = <^. 
In N-vaton scenario, usually we assume the curvature perturbation caused by inflaton is 
relatively small and can be neglected. The more general case with £ r = ^ ^ is discussed 
in Appendix A. Here we consider ( r = and Eq. ( 12.26 1 ) gives 



e 4C _ h£ fi^ e 3 e c + Q ffD -1 = 0. 
Order by order, from Eq. (12.281) we have 

JV 

c (1, = ^^,4), 



(2.28) 



i=l 



and 



where 



^(2) 



1 



4-n 



CT,D 



N 



|E^^(1 + ^)(C«) 2 -(8 + ^,d)(C (1)XS 



.4 



4 - fi CT ,D 
The total curvature perturbation up to second order is 

c = c (1) + k (2) = ^E^^e 

z t=l 



(2.29) 

(2.30) 
(2.31) 



oaN a ( N \ 2 

+ X^^^ (1 + ^(^i) -(i+a^^^te "^-] • (2 ' 32) 

Assume that the two different curvatons are uncorrelated with each other and then 



<CS(kiKS,(k 2 )> = WV^Jh)5^(k 1 + k 2 ). 
Using Eq. (12.291) . we can easily calculate the primordial power spectrum: 

N 
i=\ 

For convenience, we introduce a new parameter ctj as follows 
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(2.33) 



(2.34) 



(2.35) 



The constraint on the coefficients a» is 

A' 



E^,d«, = i. 



(2.36) 



i=l 



Similarly we can also calculate the primordial bispectrum: 

3A 3 N 



(c(k 1 )c(k 2 )c(k 3 )> = — e ^^^^^^(C(ki)cS(k 2 )(cS ) *cS ) j(k 3 )) 

i,j,k=l 

A N 

- (l + -fi CT , D )A 4 J2 ^D^D^D^n 
i,j,k,l=l 

x(e(k 1 )^(k 2 )(cS ) *O(k 3 )) 

+ 2 permutations of {ki, k 2 , k 3 }, 



where * denotes a convolution as follows 



(CS,*CS)(k) 



1 



(27T) 



d 3 qC£(q)CSl(k-q) 



'IjO \ x/ J(J j,o 



(2.37) 



(2.38) 



After straightforward calculations, we get 

(C(k 1 )C(k 2 )C(k 3 )) 



3 ^ 

— Y / nl iD ^(i + h i )-(2 + An (T , D ) 

< (2tt) 3 [P c (/c 1 )P c (A; 2 ) + 2 perms}5 3 (k 1 + k 2 + k 3 ). (2.39) 



Using the definition of f NL in Eq. (12.91) . we have 






(2.40) 



For single curvaton, the solution of Eq. (l2.36j) is a = l/Vl 2 aD and then f^i 9e = -A- (1 + 

/i) — | ^P-, where /d = Af2 CTi £). It is the same as the result in the literatures. 

To compare with the cosmological observations, we introduce a "dimensionless" power 
spectrum P^ which is defined by 



2vr 2 

(C(k 1 )C(k a )> = 4-^ 8 (ki--k a ). 



k\ 



The power spectrum of 8 ai t is given by 

P <W = 

According to Eq. fl2.17p and (12.191) . we have 

4 



2tt) 



P = —r? P = 



n 2 TJ 2 



(2.41) 



(2.42) 



(2.43) 



where 

% = dig*. (2-44) 

The value of a,- takes the form 



S"i,ol C 9vr 2 p c 



a 4 = ^P, /P c = — V s "- (2-45) 



Based on Eg. (12.341) . the amplitude of the primordial power spectrum P^ becomes 

A 2 N 
Pc = ^Y.f^Hl (2.46) 

In [3] WMAP normalization of the primordial power spectrum is 

P(,wmap = 2.457+°;°93 x iO" 9 - (2.47) 

On the other hand, the amplitude of primordial power spectrum generated by inflaton is 

which should be much smaller than P^^wmap, namely 

H* < 4.4 x l(T 4 v^M p . (2.49) 

Gravitational wave perturbation (tensor perturbation) is also generated during infla- 
tion in N-vaton model. The tensor perturbation only depends on the inflation scale and 

its amplitude is given by 

tt2 lyrl 

Pt = -±]f. (2.50) 

Usually we define a new parameter named the tensor-scalar ratio r to measure the am- 
plitude of the tensor perturbations: 

r = P T /P C . (2.51) 

So the Hubble parameter during inflation is related to the tensor-scalar ratio by 



n 



>V2 A 



H * = V2 P ( P ~ (2,52) 

Using WMAP normalization (j2.47p . we get H* = 10 _4 r2M p . If the density perturbation 
is dominated by inflaton fluctuation, we have r = 16e. In curvaton/N-vaton scenario, 
the density perturbation caused by inflaton is subdominant, and thus the inflation scale 
should be relatively low, i.e. r < 16e. 

Note: Going beyond the sudden-decay approximation, the single-curvaton model was 
studied in [10,13] in detail. We hope that one can do it for N-vaton model in the future. 
In this paper we always adopt the sudden-decay approximation. 
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3 Bound on the non-Gaussianity parameter f^L 

In this section we consider the case in which the Hubble parameter is roughly a constant 
during inflation. If e = —H/H 2 is large, the variation of inflaton is larger than the 
Planck scale [32]. Usually this kind of inflation model cannot be embedded into string 
theory [33-38]. For simplicity we also assume the values of curvatons don't evolve between 
Hubble exit during inflation and the beginning of their oscillations. So we have <Tj j0 = 
gi(ci,*) = Cj,*, and thus qi = I/ctj,* and hi = 0. Here we are interested in the case of large 
non-Gaussianity. From Eg. (12.401) . a large non-Gaussianity can be obtained if «j ^> 1, 
but fl a ,D is not necessarily required to be much smaller than 1. For example, if one or 
more coefficients «j are large enough and Q 2 . D (%i takes a finite value for fi ai ,D "C 1, /a/x 
can be large even when VL,j,d = 1 because of /jvl ~ (^. D a i) 2 /^^ it D- This case can be 
possibly achieved only for multiple curvatons. In the single-curvaton model, Q 2 D a = 1 
and /tvl ~ l/VL a £>. However whether the above conditions can be naturally realized in 
a concrete N-vaton model is still an open question and we will return to this problem in 
some future work. Here we only give a brief discussion for the case with two curvatons in 
the Appendix B. 

From now on, we only focus on the case with VL a ^ <C 1 for simplicity. In this case a 

large non-Gaussianity is also expected. Now A = 3/4, the amplitude of primordial power 

spectrum and the non-Gaussianity parameter in Eg. (12. 46 p and (I2.40p are respectively 

simplified to be 

1 N H 2 

Pc-J^^Lo^, (3.1) 

and 

5 N 

fNL = \Y,^D*l (3.2) 

Since aj is only constrained by Eg. (I2.36J) . usually we need more information if we want 
to constrain the non-Gaussianity parameter f^L- 

3.1 Lower bound on Jnl 

Let's introduce a very useful inequality 

N N / N \ 2 

i=l 1=1 \i=l / 



where ui > and Vi > for i = 1, 2, ..., iV. The equality in Eg. (I3.3[) is satisfied only when 
Ui/uj = Vi/vj for i,j = 1, 2, ..., N. Using this inequality, we immediately find 

N N / N \ 2 

Taking Eq. (12.36j) into account, we find the non-Gaussianity parameter Jnl in Eq. (13.21) is 
bounded from below, namely 

fsL > ^—- (3-5) 

The equality is satisfied when aiQ ffit D — @ which is a constant. We can easily check it. In 
this special case the solution of Eq. (l2.36j) is given by 

= 1/^,d, (3.6) 

and then 

Instituting this solution into Eq. ( 13.21) . we get 

Jnl = ^—- (3.8) 

Keeping fl a ,D fixed, the non-Gaussianity parameter f^L i n N-vaton model is not less than 
that in single-curvaton model. 

3.2 Upper bound on f NL 

In this subsection we take more information into account. Because we only focus on 
the limit of Q^d ^ 1> the radiation produced by inflaton is always dominant before 
the curvaton decay. After that curvatons oscillate around their minimums a"j = and 
their energy density decreases as a~ 3 . Once the Hubble parameter drops below r CT , the 
curvatons energy is converted into radiations. Similar to the arguments in [8, 10], the 
energy density parameter £l ai ,D at the time of curvatons decay is given by 

of /m,\2 



a ^=smM) ■ (3 - 9) 

Instituting the above equation into Eq. fl3.ll) . the amplitude of primordial power spectrum 

becomes 

H 2 N 

^ rn'r- S^ 1 (310) 
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or equivalently, 

E-^ = (24vr) 2 P c ^. (3.11) 

«=i * 

The WMAP normalization gives a constraint on Y,iLi m i°'i*- 

In this section A = 3/4, ^((Tj,*) = er^,* and #j = 1/oi,*. Eq. (j2.45p is simplified to be 

«« == T^-^' (3-12) 

and then 

CKjfL-. n = 4/ — -. (3.13) 

i - D i92yr ff v ; 

If vtii — m for z = 1,2, ...,iV, aiQ ffiD is a constant and the inequality in Eg. (13.51) is 
saturated. Now we have 

= — #, (3.14) 

192V r ' K J 



1 


H 2 Jal 


16vr 2 


P< 


! 


,D 


r 
192 V 





and 



/« - ^#- (3.15) 



In general different curvatons <7; has different mass m^. Using Eq. fl3.2p . (I3.9P and (I3.12p . 
we find the non-Gaussianity parameter /jvx takes the form 



#* A 2 



M = 3 x 10" 7 P C " 2 -—^E<<- (3-16) 

When mi = m for i = 1,2, ...,N, we can easily check that this results is the same as 

©ISD. 

How to determine the value of 0^* is a crucial problem in curvaton/N-vaton model. In 
the literatures a^ are taken as free parameters. In classical level it is correct. However 
for a scalar field \ i n de Sitter space, if its mass is much smaller than H*, its Compton 
wavelength is large compared to the curvature radius of the background H~ l and the 
gravitational effects may play a crucial role on its behavior. In [39-41] the authors ex- 
plicitly showed that the quantum fluctuation of a light scalar field x with mass m x in de 
Sitter space gives it a non-zero expectation value of \ 2 



ro (3-17) 



57r 2 m 2 



This result is reliable for a light scaler field with m x <C \/2H* in a long-lived, quasi-de 
Sitter inflation. Here we also ignore the possible corrections from the cubic, or higher- 
power terms in the curvaton potential. So the typical or average energy density of the 
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o rr4 

scalar field y is t^-%. Since the masses of curvatons are assumed to be much smaller than 

^ 167H 

iJ*, the total energy density of curvatons can be estimated as 167r2 * , which implies 

E"»?^ = -p-r- (3 - 18) 

Using the inequality fl3.3|) . Eg. (13. lip and (I3.18|) . we have 

N 3 
£m, ? *i < 6V6(iVPcr ff ) 5 ^M p 2 . (3.19) 

8=1 

We see that the non-Gaussianity parameter Jnl in Eg. ( 13.161) is bounded from above: 

f NL < 4.41 x 10" 6 P c " f iV^-^. (3.20) 

The ineguality (13.201) is saturated when these curvatons fields have the same mass: mi = 
rri2 = ... = m^v = Tt- Now aiVL au £> is a constant and the lower bound (J3.5P is also 
saturated. One point we want to stress is that Q, a> D is not kept fixed. Keeping the 
inflation scale H* (or tensor-scalar ratio r), the number of the curvatons N and the 
curvaton decay rate IV fixed, the non-Gaussianity Jnl is maximized in the case where all 
of the curvatons have the same mass. We discuss this special case in Sec. 3.4 in detail. 

3.3 Adiabatic condition 

In [21] the author pointed out that the curvaton model is free from the constraint of 
isocurvature perturbation in WMAP [3] if the cold dark matter (CDM) is not the direct 
decay product of the curvatons and CDM is generated after the curvatons decay com- 
pletely. So does N-vaton. Denote H cdm as the Hubble parameter when CDM is generated 
and thus H cdm < IV The Hubble parameter H cdm is related to the temperature T cdm at 
the epoch of CDM creation by H cdm = T^ dm /M p . Therefore 

r CT > ^p. (3.21) 

Combining with Eg. (13.201) . we find 

T cdm < 1.87 x lO^ivM^I GeV, (3.22) 



where we use Eg. (12.521) and WMAP normalization Pq = P(,wmap- In [42] the relationship 
between T cdm and the mass of CDM M cdm is roughly given by 



Lcdm — ^JJ-cdm- \6.Z6) 
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So the mass of CDM is bounded from above 

M cdm < 3.7 x 10 13 ivM f~l GeV. (3.24) 

For example, N ~ 10 3 , r ~ 10~ 4 and /jvx ~ 50, the mass of CDM is less than 3 x 10 8 
GeV. On the other hand, f^L is bounded by 1/M^ dm from above. 

3.4 The case with rm = m for i = 1,2, ..., N 

In this case the constraint coming from the amplitude of power spectrum (13. lip and the 
estimation of the total energy density of curvaton during inflation (I3.18P are respectively 
simplified to be 

4-E-l = (24vr) 2 P c ^, (3.25) 

and 

87r^m 2 
According to the above two equations, we find that curvaton mass m is related to curvaton 
decay rate r CT by 

m = 6.68 x 10' 6 p- l N—- 1-. (3.27) 

Keeping IV fixed, the mass of curvatons m in N-vaton is N times of that in single-curvaton 
model. 

In Sec. 3.2, we argue that the upper bound on the non-Gaussianity parameter /jvl in 
Eg. (13.201) is saturated when the curvatons have the same mass and now we have 

f NL = 4.41 x lU-s P p N \J^-. (3.28) 

Keeping r CT and H* fixed, larger the number of curvatons, larger the non-Gaussianity 
parameter Jnl- On the other hand, if iV is fixed, smaller T a , larger Jnl- However, similar 
to the argument in [29], the curvaton decay rate is larger than the gravitational strength 
decay rate, i.e. 

where c is supposed to be an order one coefficient which we have not known exactly. The 
curvaton decay rate cannot be arbitrary small. Substituting Eq. (l3.27p into (I3.29p . we find 
that r CT is bounded from below by the number of curvatons 

r 9/2 



Ml 



T a > 1.3 x lO-V^, W^— . (3.30) 
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The lower bound on the curvatons decay rate rises as the number of curvatons increases. 
Combing with Eq. ( 13.281) . we find Jnl is bounded by the tensor-scalar ratio from above 

f NL < 0.034 . c .p-*JV-i^ = 10 3 - C •(£)*. (3.31) 

For N — 1, our result is the same as that in [17] where we go beyond sudden-decay 
approximation and ignore the coefficient 3/87T 2 when we estimated the expectation value 
of square of the curvaton field. Here we introduce an order one coefficient c to encode the 
uncertain coefficient in the calculations. On the other hand, using Eq. (13.271) . (13.281) and 
(I3.3ip . we obtain 

c 2/3 / m/N \ ' _ 2 / m/N \ * 
fNL< T^{l.3xl0m p ) - C3 \2x 10 4 GeV J " (3 ' 32) 

Requiring the mass of each curvatons is smaller than the Hubble parameter H* leads to 
another bound on Jnl, i-e. 

f NL < 2.3 x 10 3 • c§ • ^. (3.33) 

If r > 2 x 10 4 / (c 4 iV), the constraint in Eq. 03.33]) is more restricted than that in ( I3.3ip . 
To summarize, the bound on the non-Gaussianity parameter f NL is given by 



f NL < min 



A „ 1 o q ^ i no 



10 a -c- — , 2.3 x lO" 5 -^ 



NJ ' 7VV3 



(3.34) 



We see the the constraint on /jvx i n N-vaton model is more stringent than that in single- 
curvaton model. The reason is that the energy density of each curvaton during inflation 
is roughly if 4 and thus the total energy density of curvatons in N-vaton model is much 
larger than single curvaton energy density. Larger the number of curvatons, larger Vl a ^. 
Since J'nl ~ 1/^o-,d, the non-Gaussianity is suppressed in N-vaton model by the number 
of curvatons N. A reasonable estimation of the number of curvatons in string theory 
might be 10 3 and r < 10~ 3 if we require the variation of inflaton be smaller than Planck 
scale [34]. If so, /jvx < 32 ■ c. Usually j'nl in N-vaton model should be less than 10 2 . 
For f NL > 10, r > 10~ 8 iV and m > 10 7 N GeV. Typically we have N ~ 10 3 and then 
m > 10 10 GeV, r > 10~ 5 which implies H* > 10 12 GeV. 

3.5 N-vaton vs. single-curvaton model 

According to previous discussions, the non-Gaussianity parameter f NL in N-vaton model 
is larger than that in single-curvaton model if the decay rate of curvaton is kept fixed, 
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and the maximum value of Jnl in N-vaton is obtained when all of curvatons have the 
same mass. The maximum value of Jnl is \/~N times of that in single-curvaton model. 
However now the curvaton mass in N-vaton is iV times of that in single-curvaton model. 
The requirement that the curvaton decay rate be larger than the gravitational strength 
decay rate in N-vaton model becomes much more stringent than that in single-curvaton 
model. That is why the upper bound on f^i in Eg. (13.311) is suppressed by a factor 1/7VX 
On the other hand, we consider the mass of different curvaton is quite different 
from each other. For simplicity, we estimate a^* ~ H^/rrii. According to Eg. (13. 101) 
and (I3.16p . the contributions to the amplitude of primordial power spectrum and non- 
Gaussianity parameter from curvaton ctj are respectively Pq ~ H% / '{M^Y a mi) and /jv£,i ~ 
H*f \PfMpYamf) . If the lightest curvaton a^ is much lighter than other curvatons, the 
total primordial power spectrum and non-Gaussianity are roughly contributed by a^. Now 
N-vaton model is reduced to single-curvaton model and Jnl ~ ^l/(-P<* 2 H*). Similarly, 
requiring IV > m\j{c t M'^) yields /nl < 10 3 ■ c ■ r*. 



4 Spectral index and non-Gaussianity 

The spectral index of the primordial power spectrum generated by curvatons is defined 

as 

d1nP?° 

2e + 2 V(T(n (4.1) 

ml 



n't = 


1 + 


dink 


= 1 


TV 

E^ 

i=l 


i: n a i 


1 d'< 


l V{°i) 


3H* 


daf 



where 



If the primordial power spectrum is dominated by the curvature perturbation generated 
by curvatons, n s = n™ c . The masses of curvatons are assumed to be much smaller than H* 
and then n s ~ 1— 2e. For n s = 0.96, e = 0.02 which might be realized in landscape inflation 
[43-45] or the monodromies [46]. However in this case the Hubble parameter H during 
inflation cannot be taken as a constant any more. In [24], we showed that the values of the 
curvatons depend on the initial condition of inflation which should be fine-tuned to achieve 
the suitable amplitude of primordial power spectrum and non-Gaussianity parameter Jnl- 
It is quite unnatural. Usually a closely scale-invariant curvature perturbation generated 
by curvaton is expected in curvaton/N-vaton scenario. 

On the other hand, the spectral index of the primordial power spectrum generated by 

inflaton is 

d In P inf 
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In some inflation models, e ~ 0, but the order of magnitude of rj can be — (9(10 _1 ) 
to — O(10~ 2 ). If the inflaton fluctuation makes a significant contribution to the total 
primordial power spectrum, a red-tilted primordial power spectrum is possibly obtained. 
We calculate the curvature perturbation for this scenario in Appendix A. Introduce a 
parameter (3 to measure the relative amplitude of power spectrum generated by curvatons: 



(3 = P c nc /P c to4 , (4.4) 



and then 



yinf _ / 1 n\ r>tot 



p™* = (l - (3)Pl ot . (4.5) 

Now the spectral index becomes 

dlnP;°< . f 

n s = l+ -^ = /K c + (1 - P)n7 f 

= 1 - (6 - A(3)e + 2(3 Va(7 + 2(1 - (3)r), (4.6) 

where a, in Eq. (14.21) should be replaced by 7$. We consider e<l and rj acT <C 1 and thus 

n s ~ 1 + 2(1-/3)77. (4.7) 

For (3 = 0.8 and 77 = -0.1, n s ~ 0.96. 

Since the total primordial power spectrum is not only generated by curvatons in this 
scenario, some formulations in Sec. 3 should be modified by some powers of (3: f^L is 
replaced by /nl/P 2 and WMAP normalization becomes P^ = /3P(,wmap- For example, 
Eq. (13.341) is changed to be 



f NL < min 



r \i „ 4 2 r 1 / 6 



10 d -/?4 -C- — , 2.3 x 10 d -/?3 • C 3 



N ' ^ iVV3 



(4.? 



For /3 = 0.8, the bound on the non-Gaussianity does not change so much and a large 
value of /jvx is still achieved naturally. But we need to stress that a large value of Jnl is 
obtained only when r is not too small. 

The size of non-Gaussianity generated by inflaton is controlled by a factor (1 — j3) 2 
in Eq. (lA.18p and rich phenomena are expected in this mixed scenario. A large value 
of f™[ is possibly detectable if (1 — (3) is not so small. Another bonus of this mixed 
scenario is that if the adiabatic fluctuation generated by inflaton is big compared to that 
from the curvaton {(3 ~ 0), our model is relaxed from the constraint on the isocurvature 
perturbation in [3] even in the case where dark matter was generated before the decay of 
the curvaton. 
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5 Random matrix and typical mass spectrum in string 
theory 

Axions are typically present in large numbers in string compactifications, and even when 
all other moduli are stabilized, the axion potentials remain rather flat as a consequence 
of well-known nonrenormalization theorems [47] . Following [48] the potential of N axions 
'-Pi is 

N r / 



v^) = J2Ai 



i=l 



COS . 

V Ji 



(5.1) 



where /, is the axion decay constant and A, is the dynamically generated scale of the axion 
potential that typically arises from an instanton expansion. Redefine the axion field as 
Oi = f/fi and then the Lagrangian for small axion displacements Oi <C M p in [49] is given 
by 



£ = E 



N r i i 

-{da) -^m i o i 



i=\ 



(5.2) 



In [48,49] the axion fields are taken as inflatons and the value of <Tj is larger than M p /y N, 
but smaller than M p . In [35] we argued that the vacuum expectation value of <Tj is bounded 
by M p /yN from above and this inflation model might be inconsistent with full quantum 
theory of gravity. In this section we suggest that these axion fields play the role as 
curvatons, not inflatons. 

It is still difficult to explicitly calculate the mass of axion in the context of KKLT 
moduli stabilization [50]. However, in [49] the authors found an essentially universal 
probability distribution for the mass square of axions as the Marcenko-Pastur law 



1 



J{b — m? / rn 2 ) (m 2 / 'm? — a 



P( m ) = 7^ — 5 » ( 5 - 3 ) 



2tiv m 



for a < m 2 /m 2 < b, where 



a 



1 - v^) 2 , (5.4) 



b=(l + ^) 2 . (5.5) 

The shape of the distribution only depends on a single parameter v which is determined 
by the dimensions of the Kahler and complex structure moduli spaces. This distribution 
is universal because it does not depend on specific details of the compactification, such as 
the intersection numbers, the choice of fluxes, or the location in moduli space. It is also 
insensitive to superpotential corrections. But we cannot determine the overall mass scale 
from string theory. In a KKLT compactification of Type IIB string theory, there are h\j. 

17 



axions, and hi t i + ft.2,1 + 1 is the total dimension of the moduli space (Kahler, complex 
structure, and dilaton), so that 



hi,i + h 2 ,i + 1 

In [49] the authors argued that the models of v — | are strongly favored. 

In general, the number of axions is roughly O(10 2 ~ 10 3 ) in string theory. So 

1 N rbfh 2 



(5.6) 



1 N -bjn 2 

— J2 m i k = ( m2k }— m 2k p{m 2 )dm 2 = rh 2k s(k,v) (5.7) 

iV ■ i Jam 2 



up to the order of 1/N which can be safely neglected in our analysis, where k is just a 

number and 

1 r b i 

s(k,v) = - — / x k ~ l J(b-x)(x - a)dx. (5.8) 



2TTV Ja 

The function s(k, v) has some interesting properties: 

s(0,v) = s(l,v) = s(k,0) = l. (5.9) 

Since s(l,v) = 1, (m 2 ) = fh 2 which denotes the overall mass scale. We also define 

lf>faN(- 2 V> = (- 2 V>>. (5-10) 



N l=l 



Here we have (a 2 ) = „ 2 % and then 

\ ' 8ir z m^ 

EA 2 = ^ 2(M 4-M). (5.H) 

Because s(0,v) = 1, Y^iTTifof = g7r2 * and Eq. (l3.18l) is automatically satisfied. Here 
we assume that cr, <C /j. Otherwise, the quartic, or higher-power correction terms from 
the expansion of the axion potential will be important. In string theory, the axion decay 
constant fi is generically large [51] and our assumption of a"j <C fa is reasonable. 

According to Eq. (l3.10l) and (13.161) . we can easily calculate the amplitude of primordial 
power spectrum and the non-Gaussianity parameter generated by curvatons: 

P? = 6.68xlO~ e s(-l/2,v) J ^-, (5.12) 



NH 

c 



f n N c L = 1.14 x KT 8 s(-l/4,?;)(P"r 2 ^-. (5.13) 



M$TZrh* 
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Here are three unknown scales, if*, r CT and m, which still cannot be determined by 
microscopic physics. Canceling IV, we have 



i m 



fW L = 0.66/ 1 (^(P c " c )-^iV"5 

-n* 



where 



f 1 (v)=s(-l/4,v)/sH-l/2,v). 
On the other hand, canceling fh yields 



f n N l = 4.41 x lO' b f 2 (v)(P^)-N2 



i Hi 



M*r a - 



with 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



f 2 (v)=s(-l/4,v)/s*(-l/2,v). 

We have /i(l/2) = 0.76 and / 2 (l/2) = 0.98. The behaviors of f\{v) and /2(f) are showed 
in Fig. 1. When v —* 0, the mass gap of curvatons disappears and we can expect that our 
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Figure 1: The function fi(v) and /2(f)- 

model is reduced to to the case in Sec. 3.4 where all of curvatons have the same mass. 
Since s(k, 0) = 1 and then /i(0) = ^2(0) = 1, we see both the amplitude of primordial 
power spectrum and the non-Gaussianity parameter are really the same as those in Sec. 
3.4. On the other hand, in Sec. 3.2, we find that the non-Gaussianity parameter /jvx is 
maximized when mi = m 2 = ... = m^ = m for keeping H*, T a and iV fixed. This model 
is really consistent with our analysis: /2(f) approaches its maximum value when v — > 0. 
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Here we also know how the mass scale m varies with v. For a given f^f L , m ~ l/fi(v) 
rises as v increases. 

In this case, we can also calculate rj a(T , i.e. 

1 



Va 



]t{v) 



— 2 

m 

w 



where 



t(v) = s(1/2,v)/s(-1/2,v). 
The function t(v) is illustrated in Fig. 2 and r(l/2) = 0.73. 



r(v) 




(5.1i 



(5.19) 



Figure 2: The function t(v). 

In general, curvatons fluctuations only contribute to a part of the total primordial 
power spectrum, P" c = /3P^ ot , and then f NL ~ /3 2 f^f L . Similarly, requiring T a > 
rh 3 /(c 4 M%) yields 



f NL < min 
where 



10 3 ■ /?t -c-d^v)- {jzY, 2.3 x 10 3 • /?t . c § .d 2 ( v ) 



d 1 (v) = s(-l/4,T;)/a?(-l/2,v), 

d 2 (v) = s(-l/4,v)/sl{-l/2,v). 



r l/6 

ATV3 



(5.20) 



(5.21) 
(5.22) 



The functions d\(v) and d 2 {v) are shown in Fig. 3, and di(l/2) = 0.81 and d 2 (v) = 0.79. 
Again we see that our results are exactly reduced to the case in Sec. 3.4 when the 
parameter v approaches to zero. 
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di ,( v ) 




Figure 3: The function di(v) and d 2 {v). 



5.1 Compare to experiments 

In this subsection we focus on how to compare our model to experiments. We consider 
the case with Pp c = f3Pf ot . There are 8 parameters: 



inflation : if*, e, 77 
N-vaton : N, T a , fn, v 
ratio parameter : /3 



Using Eq. (15.181) . we have -^- = J?>rj aa /T{v) and then 



f NL = 1.14/32/ 1 ( v )(P c to *)-2Ar- v / r/(T(T / r (,;). 



(5.23) 



Since P 1 / = * Sm p = (1 — /?)P/ oi , the relationship between tensor-scalar ration r and e 



is given by 

r = 16(1 -0)e. 

The spectral index is given in Eq. (14.61) as 

n s = 1 - (6 - Af3)e + 2(3r] aa + 2(1 - /3)rj. 



(5.24) 



(5.25) 



To summarize, there are four quantities which can be measured by experiments: P/ ot , 
fNL((3,v,N,r) aa ), r(P,e) and n s (P,e,r),r) a(T ). For a given inflation model, which means 
e and rj are given, the parameters 770-0-, N and /3 can be determined by experiments for 
the preferable model with v = 1/2. Furthermore, if the number of curvatons is given by 
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the string theory, we can check whether our model is consistent with experiments. For 
example, Let's consider an inflation model with e = 10~ 4 and rj = —0.1, and we also 
assume fh/H* <C 1 (or r/ CTCT C 1). The tensor perturbation is too small to be detected. 
For n s = 0.96, P^ = 2.457 x 10" 9 , f NL = 30, v = 1/2 and N = 10 3 , we find (3 = 0.8 and 
r]aa = 0.0042. Now r = 3.2 x 10" 4 , H* = 4.36 x 10 12 GeV and then fh = 5.7 x 10 11 GeV. 

6 Discussions 

In this paper we explicitly calculate the primordial curvature perturbation in N-vaton 
model. Multiplicity of light scalar fields is generic in the theories going beyond standard 
model. Even though the total energy density of these light scalar fields is subdominant 
during inflation, the perturbation produced by them can dominate the density pertur- 
bation on large scale. We also suggest a realistic N-vaton model in which the axions in 
the KKLT compactification of Type IIB string theory are taken as curvatons, and a rich 
phenomenology is shown up. If a large local-type non-Gaussianity is confirmed by the 
forthcoming experiments, it can shed a light on these light scalar fields. 

In order to fit the spectral index from WMAP data, the inflaton fluctuation is still 
required to play a significant role in the total primordial power spectrum. Generally the 
tensor-scalar ratio is required to be not smaller than 10~ 5 if f^f? 1 > 10. Many inflation 
models constructed in string theory, such as brane inflation [52,53], happen in a quite 
low energy scale with r ~ 10~ 10 which is too small to generate a large non-Gaussianity in 
curvaton/N-vaton scenario. How to construct an inflation model with r ~ 0(1O~ 5 — 10 -3 ) 
and rj ~ —0.1 is still an open question. 

In general, the curvaton decay rate mediated by particles of mass Mx is expected to 
be of order m 3 /M x . So it is natural to assume that the different curvatons have different 
decay rates and decay at different time, in particular for the case where they have different 
masses. The authors in [15] gave a concrete example to show that a large non-Gaussianity 
can be obtained even when the curvatons dominate the total energy density at the time 
of decays in the case of two curvatons. However if there are hundreds or thousands of 
curvatons, whether this enhancement of the non-Gaussianity is generic or not is still an 
open question. It is worth studying this problem in the future. 

Acknowledgments 

We would like to thank Daniel Chung, Hironobu Kihara, Han-Tao Lu, K. P. Yogendran 
and Yu-Feng Zhou for useful discussions. 



22 



A Appendix: £ r = Q ^ 

This is the most general case. We also expand the curvature perturbation ^ to second 
order as follows 

C = C?> + j<f. (A.i) 

Now Eq. (12.261) becomes 

e 4C = fe^e 3 ^) e c + (1 - fi CT , D )e 4 ^. (A.2) 

Order by order, from the above equation we have 

A' 



C {1) =A^^nei + BC { ;\ (A.3) 



and 



where 



qan , / jv \ 2 

c (2) = yE^(i + ^((S) -(2 + ^, D )^ 2 E^, D e 

z i=i \i=i / 

N „ 

- 8A 2 B£n^cg,C, w + 5 2 c (Ci 1} ) + *cf , (A.4) 



A = ^— , g= ^^ , and C=J5f£-A (A.5) 

The total curvature perturbation is given by 

C = C (1) + \C {2) - (A.6) 

We assume all of these fields including curvatons and inflaton are independent. Thus 
the two different curvatons are uncorrelated with each other 

(C£(ki)Cg,l(k 2 )> = (2nfV^ »M 3 (ki + k 2 ), (A.7) 

and curvatons are also decoupled to inflaton <fi 

(CiIi(ki)Ci 1) (k 2 )) = 0. (A.8) 

The primordial power spectrum Vu generated by inflaton is defined as 

(C ) (ki)C ) (k2)) = (2 7 r) 3 P c / 3 (k 1 + k 2 ). (A.9) 

The fluctuations of curvatons and inflaton contribute to the total primordial power spec- 
trum which is given by 

V tot = v nc + v inf^ ( A10 ^ 
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where 

N 



Vr = A*Y,^ D V^. (A.11) 



i=l 



is the total curvature perturbation generated by curvatons and 

Vf f = B 2 V (4i (A. 12) 

is the curvature perturbation generated by inflaton. For convenience, we introduce pa- 
rameters (3 and 7» as follows 

= V r l c lVf\ (A. 13) 

V Gio = A~ 2 lt Vr- (A. 14) 

Thus we have 

N 
i=\ 

When (3 = 1, all of the primordial power spectrum is generated by curvatons and 7, = ccj. 
Now we have 

V^=A-*fh(P?, (A.16) 

and 

V^ = B- 2 {\-p)Vf t . (A.17) 

Similarly, we can also calculate the total non-Gaussianity parameter f^ L . Here we only 
give the result: 



etot _ n2 f nc , a/i _ o\ fcross j_ ft _ q\2 finf 



where 



m = P'm + p(i - flfizr + (i - pyfsL (A.18) 



/^ = ^EC* + fc) - (3 + V^' d) ' (a - 19) 

20 

/yT = -yA (A.20) 

/$ = ^ + ±/£/, (A.21) 



and f™[ is determined by concrete inflation models [54-63]. 
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B Appendix: Another way to get a large non-Gaussianity 
in N-vaton model 

In this section, we consider that there are two curvatons whose masses are Vfi\ and iri2 
respectively. Without loss of the generality, we assume m\ > m 2 . Once the Hubble 
parameter drops below mi, the curvaton <j\ starts to oscillate and its energy density goes 
like ~ a~ 3 . Similarly, when H ~ m.2, the curvaton 02 begins to oscillate. For simplicity, we 
assume that the universe is dominated by radiation before 02 starts to oscillate. The scale 
factor at the time when 0\ starts to oscillate is denoted as a = 1, and then a = Jmi/m^ 
when H ~ 7712- Since the energy density of an oscillating curvaton goes like ~ a -3 , the 
ratio of the energy density between these two curvatons at the time of their decay is 

x= — --J j-, (B.l) 

Pa 2 ,D V m 2 #2,* 

and then 

x 

^o-i ,D = Z &cr,D, Q<T 2 ,D = ~ &>a,D- (B-2) 

1 + X 1 + X 

Now we have a t = a c H 2 /a 2 ^, where a c = A 2 /(9tt 2 P^) is just a numerical coefficient. The 
constraint on a, in Eq. (l2.36|) reads 

a^ D ^=(l + ^x) . (B.3) 



According to Eq. (j2.40p . the non-Gaussianity parameter f^L becomes 

5 1 (l + .x)(l + ^x) 

fNL " Un^ (ITS • (R4) 

V V m2 

If mi = m.2, /jvl is large only when fi CT] D <C 1. On the other hand, if mi ^> rri2 and 
x > 0(1), fNL ~ 1/(^o-,d) which is large when Q a ,D <C 1. Now if x <C 1, but c<j = 
m\/m,2X > 0(1), we have 

f"L ~ T7 n z- , ,.x 2 a/— (B.5) 



which can be large even when f2 CTj D ~ 1. However i< 1 and uj > 0(1) cannot be achieved 
if the values of curvatons during inflation take the typical values of* = 3Hf/(87i 2 m 2 ), 
because x = (m2/mi) 3 ' 2 and thus u = x 2 ^ 3 <C 1 if x <C 1. Here we also want to remind 
that the assumption that these two curvatons have the same decay rate is not reasonable 
if nix/ 7712 ^ 1- 
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